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Abstract: We study cosmological perturbations for a ghost free massive gravity theory formu- 
lated with a dynamical extra metric that is needed to massive deform GR. In this formulation FRW 
background solutions fall in two branches. In the dynamics of perturbations around the first branch 
solutions, no extra degree of freedom with respect to GR is present at linearized level, likewise what 
is found in the Stuckelberg formulation of massive gravity where the extra metric is flat and non 
dynamical. In the first branch, perturbations are probably strongly coupled. On the contrary, for 
perturbations around the second branch solutions all expected degrees of freedom propagate. While 
tensor and vector perturbations of the physical metric that couples with matter follow closely the 
ones of GR, scalars develop an exponential Jeans-like instability on sub-horizon scales. On the other 
hand, around a de Sitter background there is no instability. We argue that one could get rid of the 
instabilities by introducing a mirror dark matter sector minimally coupled to only the second metric. 
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1 Introduction 

Recently, there has been a renewed interest in the search of a modified theory of gravity at large 
distances through a massive deformation of GR (see for a recent review [1]). A great deal of effort 
was devoted to extend at the nonlinear level [2] the seminal work of Fierz and Pauli (FP) [3]. The 
FP theory is defined at linearized level and is plagued by a number of diseases. In particular, the 
modification of the Newtonian potential is not continuous when the mass vanishes, giving a large 
correction (25%) to the light deflection from the sun that is experimentally excluded [4]. A possible 
way to circumvent the the discontinuity problem is to suppose that [5] the linearized approximation 
breaks down near a massive object like the sun and an improved perturbative expansion. In addition, 
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FP is problematic as an effective theory. Regarding FP as a gauge tlieory where the gauge symmetry is 
broken by a explicit mass term m, one would expect a cutoff A2 ^ mg~^ — {mMpiY^^ , however the real 
cutoff is A5 = (m**Mpj)^/^ or A3 — [ra^MpiY^^ , much lower than A2 [6]. A would-be Goldstone mode 
is responsible for the extreme UV sensitivity of the FP theory, that becomes totally unreliable in the 
absence of proper UV completion. Recently it was shown that there exists a non linear completion of 
the FP theory [7] that is free of ghosts up to the fourth order [7] , avoiding the presence of the Boulware- 
Deser instability [8] . Then the propagation of only five degrees of freedom (DoF) was generalized to 
all orders in [9]; this was shown also in the Stuckelberg language in [10]. 

Quite naturally massive gravity leads to bigravity. Indeed, any massive deformation, obtained 
by adding to the Einstein-Hilbert action a non-derivative self-coupling for the metric requires the 
introduction of an additional metric g that may be a fixed external field, or be a dynamical one. 
When g is non-dynamical we are in the framework of aether-like theories; on the other hand if it is 
dynamical we enter in the realm of bigravity [11] that was originally introduced by Isham, Salam and 
Strathdee [12]. The need for a second dynamical metric also follows from rather general grounds. 
Indeed, it was shown in [13] that in the case of non singular static spherically symmetric geometry 
with the additional property that the two metrics are diagonal in the same coordinate patch, a Killing 
horizon for g must also be a Killing horizon for 5; see [14] for a concrete example. Actually it turns 
out that the off diagonal solutions show no modification of gravity at large distance [15, 16]. Also 
cosmology calls for the bigravity formulation of massive gravity. When the second metric is static there 
is no homogeneous spatially flat FRW solution [17-19], on the contrary in the bigravity formulation 
flat FRW homogeneous solutions do exist [20-22] . See also [23] for a different approach to cosmology 
in massive gravity. In this paper we study perturbations around the FRW background solutions found 
in [20]. 

The outline of the paper is the following. After a brief introduction to the bigravity formulation 
of massive gravity in section 2, the FRW solutions are reviewed in section 3. The perturbed FRW 
geometry is introduced in section 4 and the perturbed Einstein equations are given in section 5. The 
dynamics of the perturbations in various cases are studied in sections 6, 7 and 8. Section 9 contains 
our conclusions. 

2 Massive Gravity and Bigravity 

Any modification of GR that turns a massless graviton into a massive one calls for additional DoF. 
An elegant way to provide them is to work with an extra tensor g^^^. When coupled to the standard 
metric g^^, it allows to build non-trivial diff-invariant operators that lead to mass terms when expanded 
around a background. Consider the action 

j d^x[^gKMln + ^ [Ml (n-2m^V)+ L,natt] } , (2.1) 

where R{gi) are the corresponding Ricci scalars and the interaction potential is a scalar function of 
the tensor AT/^ = g^'^gau- Matter is minimally coupled to g and it is described by Lmatt- The constant 
K controls the relative size of the strength of gravitational interactions in the two sectors, while m 
sets the scale of the graviton mass. The action (2.1) brings us into the realm of bigravity theories, 
whose study started in the '60 (see [11] for early references). An action of the form (2.1) can be also 
viewed as the effective theories for the low lying Kaluza-Klein modes in brane world models [11]. The 
massive deformation is encoded in the non derivative coupling between 17^1, and the extra tensor field 
5^1/. Clearly the action is invariant under diffeomorphisms, which transform the two fields in the same 



- 2 - 



way (diagonal difFs). Taking the limit n — > oo, the second metric decouples, and gets effectively frozen 
to a fixed background value so that the "relative" diffeomorphisms are effectively broken, as far as 
the first metric is concerned. Depending on the background value of g^ii, one can explore both the 
Lorentz-invariant (LI) and the Lorentz-breaking (LB) phases of massive gravity [25] [26]. When the 
second metric is dynamical this is determined by its asymptotic properties, as discussed below. In this 
case notice that g^^ is determined by its equations of motion (for any finite AIpi) so that we will be 
working always with consistent and dynamically determined backgrounds. The role played by g^n is 
very similar to the Higgs field, its dynamical part restores gauge invariance and its background value 
determines the realization of the residual symmetries. 
The modified Einstein equations can be written as^ 

E^.+Qit = ^T^i (2.2) 

KE^ + Q2i^G; (2.3) 

where we have defined Qi and Q2 as effective energy-momentum tensors induced by the interaction 
term. The only invariant tensor that can be written without derivatives out of g and g is X^^ = gi^gav 
[11]. The ghost free potential [7]^ ^ is a special scalar function of 1^' = (\/X){^ given by 

4 

F = ^a„F„, n = 0...4, (2.4) 

n=0 

where The Vn are the symmetric polynomials of Y 

1/0 = 1 Vl^T,, V2^T^-T2, Vs^T^ -3nT2+2T3, 



(2.5) 



with T„ = tr(y"). In [27] it was shown that in the bimetric formulation the potential V is ghost free. 
We have that 

Qi^ ^m^[VSi^ - {V Y)'^] (2.6) 

Q2'i ^m'q-^/' {V'Yr,, (2.7) 

where (y')ii = dV/dY;;; and q ^ det X ^ det{g) / det{g) . 

The canonical analysis [9] shows that in general 7 DoF propagate; around a Minkowski background, 
5 can be associated to a massive spin two graviton and the remaining 2 to a massless spin two graviton. 

3 FRW Solutions in Massive Gravity 

Let us review the FRW background solutions in massive gravity [20] that are of the form 
ds^ = a2(T) {-dr^ + dr^ + dn^) = gi f.^dx^'dx" 
ds^ = w2(r) [~c^{t) dr^ + dr^ + dil^] = g2 ^,^dx^' dx" . 

It is convenient to define the standard Hubble parameters for the two metrics 

H^p-^-=Ha, H^^-=H^co, e=-. (3.2) 
dT a a Lu a 

Solutions fall in two branches depending on how the Bianchi identities are realized. 



(3.1) 



-"^When not specified, indices of tensors related with g{g) are raised/lowered with g{g) 

•^A very similar potential having the same form but with X instead of X^^'^ was considered in [15]. 
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• In branch one, ^ = ^ is constant and satisfies the following algebraic equation 

/2(C) =0 with/2(0 = 603 4 02^ + 01- (3.3) 
As a consequence, the Hubble parameter H of g and the one of g, Huj coincide and 

^ = ^(p + Ai), Ai = ^h-6f (2a3C + a2)] . (3.4) 

In this branch the effect of the mass deformation is to induce an effective cosmological constant 
and 

^^i^ with A2 = -^[6C"(2a4C + a3)+a2] ■ (3.5) 
A2 47rG K ^ ' 

This is not very surprising, indeed, the constraint (3.3), in the spherically symmetric case, leads 
to a branch of solutions with no modification of gravity, being the graviton mass zero around a 
fiat background [16]. 

• In branch two, ^ is not constant and the Bianchi identities are realized in the form 

c = §^e-^, e' = (c-i)^e with 00, (3.6) 

and 

0-1/ 2 

^ = SttG p + m2 (603 ^3 + 6a2 + 3ai C + ao) . (3.7) 
The ratio ^ of the two scale parameters satisfies the equation 



rn? 



e2 / 8 04 ^ , e 03 \ ai 2 02 3 oq 

^ ( — -2a2 1 +q — -ai 1 +^ + — -2a3e "y 



^"^^ (3.8) 



3 

In the expanding universe, the ratio ^ scales rapidly with redshift z, 1 + z = airtoday) / aiT) . Indeed, 
assuming that the mass scale m is related to the present cosmological constant as M^i oc A and 
that for matter with an equation of state p — w p, p — pq (1 + z)^^'^~^^\ we have 

8^ _ On ^3(^+1) at early times. (3.9) 

This quantity gets very large already at redshift z ^ 10, much later than the radiation era (z ~ 10^). 
Thus, eq.(3.8) can be satisfied only in two regimes: for large or small values of ^ The large ^ regime 
is physically uninteresting because c < when matter has w > [20]. In the small ^ regime, cosmology 
is very similar to the standard one; for ai > we have 

i= ^,-3(^+1) ; (3.10) 



8nG K p \G P 



p + Pg = p{l + 0[ — ]) ; (3.11) 



w,ff = w + 0[—j, c={4 + 3w)+Oi^—] . (3.12) 



^In presence of "mirror" matter with energy density p minimally coupled with the second metric g, eq. (3.8) gets 
modified as p — > p — p/k and the small-large ^ regime solutions can be evaded as soon as ~ k p/p, when m — > 0. 
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Once the matter is so diluted that p is neghgible in (3.8) the system falls in the fixed point region and 
^ is almost constant and the universe enters in a late time dS phase. The analysis is identical when 
the same spatial curvature kc is introduced in (3.f)^ for both metrics. 

Let us now discuss the differences with the frozen metric approach where g is non dynamical. 
Formally, the non-dynamical limit corresponds to k — > oo. To make contact with the existing literature, 
we take g equivalent to the Minkowski flat metric implying clearly that c and u cannot be arbitrary. 
Indeed, imposing that the Riemann curvature tensor of g vanishes we get that fcc < and 

(3.13) 



v kc 

Thus, flat FRW cosmology with frozen second metric exists only with a negative non-vanishing spatial 
curvature [17]-[20]. When (3.13) holds, Bianchi identities can be realized only within branch one, 
leading to eq.(3.3). That is why FRW perturbations in the Stuckelberg formalism are stuck into 
branch one that is rather problematic. 

4 Perturbed FRW Geometry 

Let us now consider the perturbations of the FRW background (3) 

Oiiu = gifj-iy + hi^^ , g^.^ = g2iiiy + h2fj,v ■ (4.1) 

parametrized as follows 

^100 = -2^1, ft,2 00 = -2c^A2 

hl/20^ = Cl/2. - a,i?i/2 , 9'Vi/2. - 3^1/2. = 5^/l^^l/2., = S''^h^^,/2,j =0, (4.2) 

hi/2ij = i/2i] + diVi/2 j + djVi/2i + 2didjEi/2 +'2Sij F1/2 ■ 

Spatial indices are raised/lowered using the spatial flat metric. 

Under a gauge transformation generated by the metric perturbation transforms 



where 



— — \- = — + OJ Huj ■ 
(cuj) c 



(4.3) 



6hi = a ^ {Cdagi + gi ay ^^C" + 9l t^a dyC) , 

6h2f,i, = {Cdag2f,u + 92aud^C +92f,a d^C) . 
and for the corresponding components 

M2=HK" + C°', SB2=cH"-C, 6E2 = (:, SF2^n^(:°; (4.4) 
SCi/2i = CTi ^Vi/2i = CT: 5h^^i/2ij = Q; 



(4.5) 



*The spatial curvatures must bo equal for consistency [20]. 
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In the scalar sector we have 8 fields and two independent gauge transformations, as a result we can 
form 6 independent gauge invariant scalar combinations that we chose to be 



c 



£^Ei-E2, Bi^B2- c'Bi + (1 - c^) E[ , 

where S1/2 = ^1/2 + ^1/2- The following additional gauge invariant fields will be useful to write in a 
compact form the perturbed Einstein equations 

j-i - F2 - Fi + - n^) Si , j-2 = F2 - Fi + - n^) =1 , 

B2 = B2-c^B, + {l-c^)E'^, 

Ai = c{A2-Ai) + [c{H-H^)~c']Ei, 

A2 - c{A2 - Ai) + [c{n- n^) - c'] =1 . 

The fields J^i/2, Aif2 and B2 can be expressed in terms of the ones in (4.6), as it is shown in Appendix A. 
In the matter sector, we define the following gauge invariant perturbed pressure and density 

Spgi^Sp-Eip , Spgi ^ Sp -Eip . (4.8) 

For matter, together with pressure and density perturbation, there is also the perturbed 4- velocity 
that consists of a scalar part v and a vector part Szi 

•^"^ (4.9) 
5ui — a {diV — diBi + Szi + Cu) . 

The corresponding gauge invariant quantity are defined as 

Us^v + E[, 6vi^ Sz, +Cii. (4.10) 

The conservation of the matter EMT leads to a set of differential relations; for scalar matter pertur- 
bations we have 

Sp'g,^il + w)[p {k^u,~3^[) -SnSp] ; (4.11) 
< = {-iw -l)u,n- ^ - vl/i ; (4.12) 

while for vector matter perturbations 

dv'.^Sv.iZw-l)!^.. (4.13) 

In the vector sector we have 4 fields and 1 gauge transformation; thus, we can form 3 independent 
gauge invariant vector perturbations 

^l/2^ =Ci/2.-V;/2,, X^=Cl^~C2^. (4.14) 
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5 Perturbed Einstein Equations 

Let us start with the scalar sector. The leading order perturbed Einstein equations for g are 
2A$i + m {^{H - 3>'i) + a^m^ f^i^Ti - AS) = -Sna'^GSpgi ; 



+ 871- G (p + p) Us 



0; 



(a^a,- - A) (a^ /i m^f: - $1 - *i) + Sij [m^ (2 A J^i + /a ^1) + 2^-1 (U^ + 2 W) 

-2 - 2 -H (2 - = SttG (5^ Jpg; , 



where 



/i =e [2e (3a3ce + a2(c+l))+ai] , h^i (6a3C'+4a2^ + ai) 
For the metric g we have 

2 2 

2 c^A^s + 6W„ (-J-sH^ - $2) + °, (Af - 37-2) = ; 



2 0(^2^. -^-2)- . ±'^1. ^2 



0; 



-c - Jjj A) 



f + c ($2 + ^2) 



+ 4 



2 2 
TO a 



(2c/iJ-2 + /2^2) + 



2 ^2 + 2K - 2 - ^c. *2 - 2*2 + 2 - - 2W„ $'2 + 2^c. *2 



= 0. 



For the vector sector the perturbed Einstein equations are 

AFi i „ _ , . ^ . TO^ 



- SttG (p + p) - — — /2 Xi = ; 



2a2 — ■ ' (i + c) 

+ 2?^ = TO^ /i 5(iVi2,) ; 

AT/2. , to2/2 



2a2c (l + c)^^ 



2Xi = 0; 



%^2,) + 



771 Qj C f\ 

%^2i) + ^^2 ^(^^^3) = ; 



where 



Vl2i=Vu-V2i, Vi2i = Vu-V2i. 

Notice that Vi2i = Xi — Vi2i- 

Finally, for the tensor perturbations we obtain 



TT _ uTT\ _ r, . 
ij "■2ijj - 



2ij 



+ 



2 [H + 



j'rp m? /i C 2 (uTT uTT\ 



(5.1) 
(5.2) 

(5.3) 
(5.4) 

(5.5) 
(5.6) 



(5.7) 

(5.8) 
(5.9) 
(5.10) 

(5.11) 
(5.12) 

(5.13) 



h 2ij -c Ah 2ij a [hiij - h^ij) = . (5.14) 
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Now that we have at our disposal the full set of equations, we can study perturbations around back- 
ground solutions in the branch one and two. In the following we will often use the Fourier transform 
of perturbations with the respect to a;*, the corresponding 3- momentum will be fc* and = k^ki. To 
keep notation as simple as possible the symbol of the Fourier transform will be understood. 



6 Branch One Perturbations 

In this case ^ = ^ is a non-vanishing constant such that /2(C) = 0, then T-L^ — %. From (5.6) we can 
express 4*2 in terms of <I>2, then from (5.5) we get that 5*2 = $2 = 0. Now, using (5.7) we have that 
E = J-2 = Q and from the relations in Appendix A, also Ti = 0. At this point it is straightforward to 
show that using the equations for the scalar perturbations of g we get 

[-W A -I- (Sw + 1)%^ + 2H'] + 3{w + 1)-H*'i + = (6.1) 

Not surprisingly, this equation describes the very same perturbations of GR in the presence of a fluid 
with an equation of state w. 

Also for vectors, being /2 — 0, again we have the very same equations as in GR 

AVi,-16TTGa^ {p + p)Sv,^0; (6.2) 

%F/^)-h2H%yij) =0 (6.3) 

with V2i = Vi2i = 0. Clearly, as in GR no vector propagates. In the tensor sector four modes propa- 
gates. From the canonical analysis, 8 DOF are expected but only 5 = 1-1-4 are accounted for. Thus, a 
scalar plus a vector are not present and are probably strongly coupled, at least around a FRW back- 
ground. Strong coupling was also found in the Stuckelberg approach [24] . This is not very surprising, 
condition (3.3) in flat space is equivalent to set to zero the graviton mass and, as a consequence, 
both spherically symmetric and FRW branch one solutions show no gravity modification [16, 20] and 
the extra DoF are frozen. Interestingly enough, in the Stuckelberg approach where the second met- 
ric is non-dynamical, only the branch one is available and strong coupling is unavoidable, another 
manifestation of the rather constrained nature of a theory with a priori given metric. 

7 Branch Two: Perturbations in dS Phase 

Before analyzing branch two case in full generality, it is instructive to consider a particular limit of it: 
a de Sitter (dS) background, for which we have 

p = const. ^ = const. c— I , fi ~ f2, "Hu: = H = H a , (7.1) 



-^Jl2f (204^ + 03) +2a2 + 4 
/3 K V i 



The dS phase is a fixed point of the FRW geometry of the branch two solution [20] . Introducing 

(7.2) 
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after some tedious computations one can show that all the equations in the scalar sector are equivalent 
to a single second order equation for 



+ 2'H 



2k^ 



9a2 ^2 ^2 ^jt4_l8^4 



+ 



0; 



(7.3) 



where 



m; 



I = m Vi I ^ + 1 



1 



(7.4) 



is the mass of the scalar field We have also replaced all space derivatives dm with ifc™ and 

All remaining fields can be written in terms of $_: 

1 



$+ = 1'+ = 0, 



^ $_(K$2^1^(2fc2+3a2m*) 2A;2£: ^ . 
^1 = ^-^5 9 ^o^, stt;- + ^ c. ; 



3a2 nml^^n 



(7.5) 



5 = 



2 a2 K ^2 ^2 ^2 -^2 ^2 _^ /j4 _ -^4) 



As a result, just a single scalar DoF propagates. 

For what concerns the vector sector, all vectors can be expressed in terms of Vi2j 



V2i = 



ml a2 Via, 



(l+e2K)(a2m2 +fc2) ' 
with Vi2 i satisfying a second order equation 



Vn = -f=rv2i, 



fc2 + a? m\ 



(7.6) 



(7.7) 



In the tensor sector two modes are propagating (4 DoF). The combination h^^^ = hjfj + ^2 nK^Jj is 
massless and satisfies 

(7.8) 



TT 



2nh\. 



u2 i,TT 



0. 



The previous equation is the same for tensor perturbations in GR. While for the orthogonal combina- 
tfon h^j = hffj -^K hlJj, we get 



hZ- +2Hhll' + (fc2 + a2 ml) h^l + a2 ml^^^hH = 0. 



(7.9) 



Summarizing, in the dS phase we have one scalar, one vector and two tensors that propagate, for a 
total of Is + 2y + 2t + 2t = 7 DoF, showing that all the expected DoF from the canonical analysis are 
propagating at perturbative level. We stress that this is not the case for the branch one perturbations. 
As we will see in the next section, the dynamics of branch two is similar except the presence of a 
matter fluid provides an additional scalar DoF. 
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8 Branch Two Perturbations 



The dynamics of branch two perturbations is similar to the dS phase, only more involved being ^ not 
constant and c ^ 1. Due to the complexity of the equations it is difficult and physically not very 
interesting to study them for generic values of ^. As shown in [20] and summarized in section 3, from 
early times to redshift of order one, the massive gravity FRW background solutions are characterized 
by a small value of As a result, when ^ << 1, the background solutions can be expanded in series of 
the dimensionless ratio t/tu <^ 1, where tjj is the age of the universe in conformal time. For instance, 
in the radiation dominated era 



Tu 10 \Tu J 



(8.1) 



where e — g^Q^^ = | rn? ^. In Appendix B, we give the explicit expressions for the scale factor a 
including the leading and next to leading terms for both a radiation and matter dominated universe. 

An interesting aspect to discuss is the m — )• limit. Naively, taking the formal limit to — > 
in the equations of motion of section (5) we get the perturbed Einstein equations for two separated 
GR copies. Among the perturbations of metric g that is coupled with matter, one scalar (induced 
by the presence matter) plus two transverse and tracclcss tensors propagate. In the sector of the 
perturbations of g that have no matter sources, only two tensor modes propagate. However, when we 
take into account that, using the background equations, ^ has a non trivial dependence on to, the very 
same limit is less straightforward. Indeed, the branch two background solutions are characterized by 
a value of ^ proportional ^ to to^ and this completely changes the nature of the limit to — > 0. 

In the equations for the g perturbations, see for instance (5.1-5.3), the coupling with the ones of 
g is through the effective coupling: to^ fi. In the small ^ regime, we have 

2 4 

fi Ki aim^ i^i — t- . (8.2) 

Thus, in the small to limit, the perturbations of the metric g precisely coincide with the corresponding 
in GR. For the equations that govern the perturbation of g, see for instance (5.5-5.7), the effective 
coupling with g is ™ °2 ■ Then we have that 



TTT^ ~ z Kia SirGpni finite. (8.3) 

As a result, g perturbations are not GR like in the m — )■ limit, moreover, as we will show they exhibit 
a non trivial structure in the momentum k that is very different from the linear structure in of GR. 
This peculiar behaviour stems from the interplay of the branch two background and the structure of 
perturbed Einstein equations; we will return later on this point. 

8.1 Scalar perturbations 

In the scalar sector, all the fields £, Bi and ^1/2 are non dynamical and can be expressed in terms of 
$1/2 that satisfy two second order equations; thus 2 scalar DoF propagate. Let us consider the case 



^If we believe that our theory is origin of Dark Energy we have to take as reference value m ^ 10 eV ~ Hq so 
that e ~ ^ . 

® Recall that at leading order f = §^4?— . 
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of a radiation dominated universe. The equations of motion for the two propagating scalars have the 
following structure 

< + -P^b + $b = a, 6=1,2, (8.4) 

where T^ab and Mab are functions of the following dimensionless arguments , e and x = kr. Note 
that well in the radiation era (and also matter era ), we have <C 1; thus, we can expand such a 
complicated expressions obtaining at leading order 

4 / T \ 

'i'i' + - <I''i + ^ _ =0; (8.5) 

T 3 \Tu J 

10a;2+42 . -5 a;^ - 15 .x'^ + 333 + 999 ^ 

<i>o + $n H $9 — 



(a;2 + 3) ' 3t2 (a;2 + 3) 

36 3(5.2+39) (8.6) 



T (2:2 + 3) ^ r2 (a;2 + 3) ' \tu 

The full expressions are lengthy and not particularly illuminating and are given in Appendix E up to 
the next to the leading order. We note the various functions 2? and admit a formal expansion in 
power of e (i.e m) equivalent on dimensional grounds to an expansion in power of t Jtu . Clearly, the 
equation for i>i is the same than GR plus small corrections, while the equation for 3'2 has a non trivial 
kr structure whose origin is the effective coupling with the metric g proportional to rr? that does 
not vanish in the limit m — > 0. In order to get some physical insight, let us consider the case x « 1 
that physically corresponds to modes well outside the horizon. We have that 

$1 const; (8.7) 

14 37 $o - 39 $1 39 
$'2' + ^ + = => '^2^^^!= const. (8.8) 

T 37 

Thus both scalar perturbations are frozen for the modes well outside the horizon. On the other hand, 
in the opposite limit, x » 1, e.g. for the modes well inside the horizon, we have 

$i^^cosfcr; (8.9) 

$'2' + '-^^ - ^$2 - ^$'r - . , ^ $2 - (fc^^"'^'" ' ' + ^('^^^ • 

The solution for the homogeneous part of the $2 equation has runaway exponentially behavior. Such a 
tachyonic instability is due to the sign of the coefficient of $2 , positive for super-horizon perturbations 
(and then stable) and negative for sub- horizon perturbation (and then unstable). For the matter 
dominate case (see Appendix D.l) the situation is very similar and the same kind of instability for the 
sub-horizon modes is present. Notice that such an instability is not present in dS case, see eq. (7.4). 

The leading contribution in the coefficient of $2 for sub- horizon modes ( a; >> 1 ) can be computed 
for a generic equation of state w, in the small ^ limit. The result is 



M 



22| 



_ ^ 2{a,[w + l)n-2a,) ^ ^ ^ ' 

fli 



fc2. (8.11) 



Thus, sub-horizon instabilities are present only when w > — ^ and the $2 perturbation grows expo- 
nentially as e(i+2'")'='^. So, for w > — i, the exponential growth of $2 invalidate perturbation theory 
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at time T ^ 1/k. As a consequences, already in the radiation dominated era sub-horizon perturbations 
become non perturbative, in sharp contrast with GR where matter perturbations become large only 
when the universe is non relativistic due to Jeans instability. The other scalar fields are given as a 
function of $i 2 in Appendix C. 

8.2 Vector perturbations 

In the vector sector, using (4.13), as in GR, the velocity perturbation can be easily obtained 

(5f, ^J^o.a'"-', (8.12) 

with Svoi an arbitrary function of k. From (5.8-5.11) one can show that all vectors can be expressed 
in terms of V12 that satisfies a second order equation given in Appendix F. Thus, only the vector V12 
propagates. For instance, in the case of a radiation dominated universe we have at the leading order 



Sv — 6vo{k) =^ constant in time; 



V2 = 



(F r2 + 5)^12 
8fcV^ + 50 



40 



k^ r2 (P r2 
3 



5) 



6vo : 



2 (fc2 + 5) V12 - 



320 



Then for super horizon modes with fc r ^ 1 we get 



V12 = Cl T- 



C- 



2T 



i (721-9) 



fc2 



64 
77 



16tA:2 



125 



(8.13) 
(8.14) 

(8.15) 
(8.16) 

(8.17) 



where Ci_2 are arbitrary functions of k. For any reasonable choice of 5vq there is no growing mode. 
The structure of the equations is similar in the case of matter dominated universe. 

8.3 Tensor perturbations 

The evolution equations for tensor perturbations in the radiation era at next to leading order are 

'4ao 



+ k'^ hi+e 



5 T 



0; 



10 



h'^ + 25 fc2 hi 



15 



{hi 



4 (ao K, — 36 02) 

5 KT 



h'. 



h2) + 
162 02 



(8.18) 
(8.19) 



hi 



/ 162 02 40fc2 (aoK- 6a2) 
V K r2 K 



= 0. 



Tensor perturbations hi of g behave as in GR, while the one of g, /12, beside a sizable coupling with 
hi at early times, show a larger damping factor (^ — > ^) and an effective larger mass {k'^ — > 25 k^). 
In figure 1 we show the numerical solution of (8.18, 8.19) . Clearly there are two regimes, depending 
on the value of kr. At very early times, when fc r <C 1, the two fields hi 2 are almost equal and 
constant due to the large coupling proportional to {hi — ^12)/'''^ present in (8.19). As soon as A; r ^ 1, 
hi starts to oscillate and triggers the oscillations of the /12 sector (and not the opposite!). Notice that 
the damping of /12 is much larger than the one oi hi, indeed 

1 , .. , 1 



hi 



while ho 



(8.20) 
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Figure 1. Tensor perturbations hi ^2 in a radiation dominated universe. The initial conditions are set in deep 
radiation era = 0.2, H ^2 10"''^ GeV) at z = 10^° as hi = 1, /12 = 1/4, h'1^2 = 0. 

During the matter era, the results are similar. 

Summarizing, hi follows closely GR perturbations; on supcrhorizon scales {kr ^ 1) /12 is propor- 
tional to hi, while sub-horizon modes are greatly suppressed (/i2 ^ hi when kr ^ 1). No instabilities 
are present. 

9 Conclusions 

In this paper we studied cosmological perturbations around FRW background solutions in a non- 
linear ghost free massive gravity theory. To construct a massive deformation an auxiliary metric is 
mandatory. In the Stuckelberg approach the extra metric is taken to be a non-dynamical Minkowski 
metric. Besides having an "absolute" object that take us back to the time of aether-like theories, 
such a theories are rather rigid. For instance, there is no black hole solution featuring modified 
Newton potential [13, 16] and there is no spatially flat FRW solution [17]-[20], thus a spatially negative 
curvature is required [19, 20]. Even if we allow spatial curvature, the resulting solutions feature an 
effectively turned off graviton mass [16, 20]. Thus, it is not very surprising that we found the 
perturbations of the branch one are exactly the same as GR. The extra modes present in the theory 
have zero kinetic term and do not propagate in a FRW background, this is typical of strong coupling 
and their propagation is expected at higher order. Perturbation theory around branch one solutions 
cannot be trusted due to the presence of the strongly coupled extra modes. In the Stuckelberg approach 
branch one solutions are the only available and nothing more can be said, at least using perturbation 
theory. 

When the second metric is dynamical, the branch two opens up. In this case perturbations are 
much richer: all 7+1 expected modes propagate, and the reliability of perturbation theory results is 
a non-trivial issue. In de Sitter all the extra modes propagate and no instability is present at all the 
length scales. In matter or radiation dominated period, the background solutions is characterized by a 
small S, value proportional to the graviton mass m. Such a dependence is the origin of deviation from 
GR in m limit with a non trivial k dependence. In particular, the perturbations of the metric g 
that couples with matter are very similar to the corresponding ones in GR, while g perturbations are 
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not GR-like in the m ^ limit. Specifically, we found well behaved perturbations for all vector and 
tensor modes, while one scalar mode shows an exponential Jeans-like instability already well in the 

radiation epoch, as soon as its wave length enters inside the horizon. Thus, though g perturbations 
are well behaved, such a fast growth drives sub-horizon scalar perturbations into a non perturbative 
regime just after a few Hubble times. 

A possible way out might be the introduction of a mirror (dark) matter sector minimally coupled 
to the metric g. As pointed out in the footnote 3, the small ^ regime is the source of instabilities and 
it can be avoided by the presence of mirror matter reestablishing the normal power counting for small 
m. The pressure of mirror fluid can contrast the development of sub-horizon instabilities. Of course, 
the analysis of such a bi-metric theory plus matter and mirror matter is a totally different ballpark 
and it requires a dedicated investigation that will not be given here, we limit to stress that such an 
approach would probably not pass the Occam's razor test. 

Acknowledgments. We thank F. Nesti for many interesting discussions. 



A Useful Relations 

One can verify that the following relations among the additional gauge invariant scalars hold 

B2 = Si - (1 - c^)£' ; (A.l) 

H^T^-HTx = {H-H^){^i-^2); (A.2) 

c\T2 + Ti) = (Bi - e){U + U^) - 2c'($i - $2) ; (A.3) 

c\A2 - ^1) = (61 - £') [c{n -n^)-d]; (A.4) 

{A2 +Ai)^2d^{^2--^i)+Bx [c {n^ + -H) - 3 c'] + 3 c' £' (A.5) 

+c [2B[ -28" - S' {n^ + n)] . (A.6) 

B FRW Background 

• During matter era, in the small ^ regime, we have p = po/a"^, po = g^^^a and for t/tjj <C 1 



t2 



,^ 7 T° Tlf 637k 



During radiation era, in the small f regime, we have P = Po = a n ^ ^iid for t/tu <^ 1 



Q=— +e ^^_g -l-e ^ (B.2) 
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C Scalar Perturbations in the Radiation Era 



For completeness we give also the leading expressions for the remaining scalars that can be expressed 
in terms of $1/2 



2)£ 

*i = -$i, ^f2 = -— -$2; (C.l) 



t., = ^t(t$;+$i), ^ = ^$i(fc2r2+3)+2T*;; (C.2) 
B2 = -60- -20t (r$^ + 5$2) ; (C.3) 

T 

^_ r2 (-9$i (fc2r2 + 9)+$2 (2A:'*T4 + 15A;2r2 + 99)+9r ($'2-3$i)) 
~ 3 (P r2 + 3)^ ■ 

D Scalar Perturbations in the Matter Era 

Leading order of the scalar evolution equations during matter era 

$;' + ^$; = 0; (D.i) 

T 

., 4(4a:4+81a;2 + 720) . -a;^ + 18 a;^ + 960 a;^ + 11808^ 

$o H $0 H $0 - (D 2) 

2 ^ r (a;4 + 18a;2 + 144) 2+ (a;^ + 18 a;^ + 144) ^ ^ ^ 

6 (a;'' + 36 .t2 + 432) ^, + 78 + 1728 a;^ + 14688^ 
T (a;'i + 18a;2 + 144) ^ (a;" + 18 a;2 + 144) ^~ 
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E Next to Leading Corrections for Scalars in Radiation 

For completeness here we give the next to leading correction for the scalar equations of motion (8.4). 

„4 o 8 4 (2alK(x^ + 3)^ + 25al (4 a;"* + 132 + 117)) 

Dii = 4 + e^ — !^ ^ ^ 2 ^; (E.l) 

4 45 T» 6075 k {x^ + sf 

„ r« 2a? (5a;'' + 6a;2 + 27) , , 

9 (a;2 + 3)' K 



V2 



36 4 (3 ao K (x^ + 9a;2 - 72) + 10 as (2a;6 + 24a;4 + OOx^ + 8I) ) 

2^2 1 — — 7: + ^^ 5 ; (E.3) 

+ 3 15 k (a;2 + 3)^ 

10a;2 + 42 



a;2 + 3 

t4 4ao (5x6 ^_ 42^4 ^Qg^i ^ 35;^^ ^ 40 (4^6 + 362;^ + lOBx^ + 8I) 

45K(a;2 + 3)^ ' 
2 



4ao 2 4agK (x2 + 3) +45a? (3x4 + 38x2 + 15) 



3 r4 9x2 ^8 405kx2 (x2 + 3)' 

„ r« a2 (2x6 + 83x^^ + 96x2 +459) 

A^i2 = e^^^-^ 2 ^; E.6 

9 (.x2 + 3)2 K ' 



■^21 = 15+^ (E.7) 
r-* 36 ao k (x2 - 27) (x2 + 9) + 10 03 (8 x^ + 177 x^ + 1305 x^ + 3807 x2 + 3159) 
~ ^ ^ 45 k (x2+3)^ 

333 - 5x'^ T* Aao (I0x« + 75x6- 54x''- 1323x2 ^ _ 
+5 02 (-16xS + 51x6 ^ i899a.4 + 8i8ix2 + 7533) 



405 K (x2 +3)^ 



(E.lO) 



F Vector Perturbations 



Equation of motion of the vector propagating mode 



2f2Kk^fn ,, 2kA:2^2^ , 9 r r , ^ 

^ J V{'2 + + h fc' V12 +No5v = Q; (F.l) 



where 

Ni = 20^/2^2 [2-^3 (<'' + e^) + 4?^^^' + ee'^H' + -H (3e'2 - CD] 

+ [/2 (4$^'^' + 8^2^^ ^ ^^^^^ ^ ^ ^^^2 _ ^^n^ ^ ^ iUf^H^] , 

and 

j'^No = 167r Gm^ 02^ ^2 (y; + 1) |2a2/2„^2 (^ (2^^ + H') + C') 
+K /2 fc2 $2 [4^2^3 _^ + CC'-H' + % (C'2 - eC")] + nk^i^nn^f^} . 



(F.2) 



(F.3) 
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All vectors are expressed in terms of V12, indeed, omitting the spatial indices, we get 
K = 2m''a''f2 [!^' + H [k f + ^)] +Kk''i^H^; 

and 



l&-KGa'^{w + l)p5v [2m'^ f2{ciH) + eUi] 2 Km^ /2 ^ 



Kk^ K ^ (p_5) 

K Kk^ 
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